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1. Introduction 

The flavor-changing neutral current transitions b — > 
s and b — > d cannot proceed at the tree level in 
the Standard Model (SM). Radiative decays such as 
B — > A s 7 , B — > X,j7 and £? — » hence proceed 

primarily via one-loop diagrams known as penguin 
and box diagrams. Because of this, these processes 
are rare in the SM, with branching ratios on the order 
of 10 -4 to 10 -6 . However, for the same reason, these 
channels are very sensitive to new physics beyond the 
SM. 

In the past few years, the B Factory detectors, Belle 
at the KEKB collider at KEK, Japan and BABAR at 
the PEP-II collider at SLAC, U.S.A., have recorded 
a total of more than 0.9 ab -1 of data, making it pos- 
sible to study these decays extensively. Both exper- 
iments benefit from high-granularity electromagnetic 
calorimetry with excellent angular and energy reso- 
lution (~ 3% for 1 GeV photons), silicon vertex de- 
tectors providing position resolution on the order of 
100 jttm along the beam axis and K/t: separation for 
momenta up to 3-4 GeV. Details of the detectors can 
be found in references [1]. 

The relevant theoretical framework for all the de- 
cays reviewed in this paper is the operator product 
expansion (OPE) [2], with which the heavy particles 
(the gauge bosons and the top quark in the SM) are 
integrated out to obtain an effective low-energy the- 
ory with five quarks. The long-distance contributions 
are captured in operator matrix elements, whereas 
the short-distance physics is described by the effec- 
tive coupling constants, known as Wilson coefficients. 
Physics beyond the SM can manifest in new operators, 
or in vastly different values for the Wilson coefficients 
of the existing operators, compared to the SM predic- 
tions. 

The OPE framework provides parton-level predic- 
tions, so the connection to the experimental ob- 
servables require further non-perturbative corrections. 
Form factors, calculations of which introduce signifi- 
cant theoretical uncertainties, are needed to study ex- 
clusive final states (B — > K* 7 , B — > KM, etc.). On 
the other hand, for inclusive decays, heavy-quark ex- 
pansion (HQE) and quark-hadron duality keep non- 
perturbative corrections under control. This allows 
more direct interpretations of the measurements of the 
inclusive decay rates. However, exclusive decays, be- 



ing experimentally cleaner, are observed earlier and 
measured with better precision. Ratios of exclusive- 
decay observables, in which the form- factor correc- 
tions partially or completely cancel, can also provide 
complimentary information to that obtained from in- 
clusive decays. 



2. b 



S7 



2.1. Exclusive b 



S7 



Since its first observation in the exclusive mode 
B -> K* (892)7 in 1993 [3], b -> sj has been the most 
extensively studied radiative penguin process. With 
the large data samples of the B Factories, many other 
exclusive modes have since become accessible and so 
far the branching fractions of nine charged and six 
neutral final states have been measured at 3cr or bet- 
ter significance (Fig. 1) [4]. For some of these modes, 
measurements of the charge asymmetry and the time- 
dependent CP violation have also been done [5]; the 
results are so far consistent with the predictions from 
the SM. 

It is worth highligthing a few of the most recent re- 
sults. The evidence for the first baryonic final state, 
B + — ► pA'y, has been reported by the Belle collabo- 
ration in 2004 [6] . While this decay has currently the 
smallest branching fraction ((2.16±g;f|±0.20) x 10~ 6 ) 
among all the b — > S7 final states measured at 3cr or 
better significance, its penguin nature becomes appar- 
ent due to the stringent upper limit on the two-body 
decay B+ -> pA [7], 

More recently, the observation by the Belle Col- 
laboration of B + — > ^7^(1270)7, the first radia- 
tive -B-meson decay involving an axial-vector reso- 
nance, with an unexpectedly high branching fraction 
((43 ± 9 ± 9) x 10~ 6 ) has attracted some interest [8]. 
Inclusive measurements of four different Kixix^j final 
states by the BABAR Collaboration yield similarly high 
branching fractions [9]. With more statistics, the de- 
cays B — ► Kir + ir° 7 will permit the measurement of 
the photon polarization, which is almost completely 
left-handed in the SM [10]. 

Finally, the most recent development in the exclu- 
sive b — ► sj front is the confirmation by BABAR of 
Belle's earlier measurement of the B — * Knj de- 
cay [11]. Both charged and neutral modes are now es- 
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tablished at better than 5<r significance. Additionally 
BABAR has also obtained the first upper limits on the 
decay B — > K r\' 7. These upper limits (at 90% C.L.) 
and the new combined averages for the Knj final 
states by the Heavy Flavor Averaging Group (HFAG) 
are given in Eq. 1: 



2.2. Inclusive b 



B{B+ - 


-> K + nj) 


= (9.3 ±1.1) x 10~ 6 , 


B(B° 


-» #°r?7) 


= (10.3+^) x 10" 6 , 


B(B + - 


► K+n'-f) 


< 4.2 x icr 6 , 


B(B° - 


+ If ?/ 7) 


< 6.6 x icr 6 . (1 



These measurements will be useful in reducing the 
systematic uncertainties in the inclusive measure- 
ments of B — > A s 7 , and in searches for higher-mass 
K* resonances. Moreover, it is possible to measure 
time-dependent CP asymmetry in the neutral final 
states. 
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Inclusive measurement of the B — > X s j total rate, 1 
where X s refers to any hadronic system with unit 
strangeness, gives one of the most widely quoted con- 
straints on physics beyond the SM. Additionally, the 
shape of the photon energy spectrum, which is mostly 
insensitive to new physics, provides essential input for 
measuring Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix elements V c b and V u b from semileptonic B-meson 
decays. Therefore, this decay has been measured mul- 
tiple times using independent techniques. 

The current experimental world average reported by 
HFAG, (3.55 ±0.24t°;°^± 0.03) x 10" 4 , where the er- 
rors are combined statistical and systematic, system- 
atic due to the shape function, and systematic due to 
the subtraction of contribution, incorporates five 
branching fraction measurements from CLEO, Belle 
and BABAR Collaborations [4] . All five measurements 
are in very good agreement with each other and with 
the SM prediction of (3.61±g'j£) x 10~ 4 , computed at 
the next-to-leading-log precision [13]. This remark- 
able agreement implies strong constraints on various 
models of new physics. For example, in type-II two- 
Higgs-doublet model, which represents a good approx- 
imation for gauge-mediated supersymmetric models 
with large tan j3, where the charged Higgs contribu- 
tion dominates the chargino contribution, the current 
experimental value implies a minimimum charged- 
Higgs-boson mass of 350 GeV [14]. 

The moments obtained from the photon energy 
spectra are consistent among the different measure- 
ments and with the HQE predictions based on inputs 
from the hadron-mass and lepton-energy spectra mea- 
sured in B — > X c lv decays. A fit to all the avail- 
able data from B — * A s 7 and B — > X c £v decays 
has yielded precise values for the HQE parameters nib 
and fi%, as shown in Table I. Further details on the 



Table I Heavy-quark expansion parameters m/, and /x^., 
and their correlation p, obtained from the combined fit to 
B — ► Xs'yanA B — > X c ^data in two HQE schemes [15]. 



Scheme 


mi, (GeV) 


/4(GeV 2 ) 


P 


Kinetic 


4.590 ± 0.039 


0.401 ± 0.040 


-0.39 


Shape Function 


4.604 ± 0.038 


0.189 ±0.038 


-0.23 



derivation of these parameters and the extraction of 
\V c b\ and \V u b\ can be found elsewhere in the FPCP 
conference proceedings [16]. 



Figure 1: Branching fractions for b — > s/d photon pen- 
guins [4]. The value for inclusive S7 has been reduced by 

an order of magnitude to fit it in the figure. 1 Both within the theory community and among the experi- 

mentalists, the B — » X s -f total rate usually refers to the decay 
rate for photon energy above 1.6 GeV [12]. 
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3. b 



c?7 



Belle result: 



Compared to b — > sj , the amplitude for the pro- 
cess b ^ dj is suppressed by the ratio of the CKM 
matrix elements |Vtd|/|Vt s |. Therefore it is experi- 
mentally more difficult to observe than b — > 57 : the 
branching fractions for even the most prominant ex- 
clusive final states are expected to be on the order of 
10 -6 . Furthermore, misidentified b — > S7 decays add 
up to the other substantial backgrounds (mostly high- 
energy photons produced through 7r° Jt\ — ► 77 decays 
or via initial-state radiation in e + e~ — > qq continuum 
events, where q — u,d,s,c) that are common to the 
measurements of all the radiative penguin decays. 

The search for b — > d/y has been carried out by 
CLEO, BABAR and Belle collaborations in one charged 
(B + — > p +r y) and three neutral (B° — » p°/uj/(f)j) 
exclusive decays. It is customary to combine the 
charged and neutral p and u> modes using isospin 
relations to obtain an average branching fraction 
defined by: B (B — ► (p/w)"/) = \{B(B+ -> 
p+7) + 7^ - A)+B(*° - 0,7)]}, 

where is the ratio of the B-meson lifetimes. The 
branching fractions for each mode from Belle and 
BABAR are listed in Table II. For the p and u> modes, 
center values are given for comparison, even in the 
absence of any significant signal. While Belle and 



Table II Branching fraction and significance for each b 
d-y mode from Belle [17] and BABAR [18]. 



(p/u,) 7 ) = 1.32± j^° 



(2) 



Mode 


B(icr 6 ) 


Sig. 


B(1(T 6 ) 


Sig. 


B+ -> p+7 

B° -» A 

B° -> w 7 
B -> (p/w)7 


n cc+0.42+0.09 
u - JJ -0.36-0.08 
1 oc+0.37+0.07 

^■^-o.si-o.oe 

n c«+0.34+0.05 
U.OU_ .27-0.10 
1 no+0.34+0.10 
J - J *-0.31 -0.09 


1.6 
5.2 
2.3 
5.1 


0.9+° l ±0.1 
0.0 ±0.2 ±0.1 
0.5 ±0.3 ±0.1 
0.6 ±0.3 ±0.1 


1.9 
0.0 
1.5 
2.1 


B -> 07 




< 0.85 (90% C.L.) 
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2, 
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Figure 2: Projections of the fit results to Mbc and AB dis- 
tributions for the Belle measurement of the decay mode 
B° — > p°7 • Curves show the signal (dashed), con- 
tinuum (dotted), B -> if* (892)7 (dot-dashed), other-B- 
decay background (dot-dot-dashed) components, and the 
total fit result (solid) [17]. 



The UTfit Collaboration has performed a fit to the 
data from the B Factories to extract |Vtd|/|Vt s | from 
the ratio B (B —> (p/u)j)/B(B -> K*j). To rep- 
resent the theory uncertainties in the ratio of the 
QCD factorization expressions for B — > (p/lu)j and 
B —> K* 7 decays and in the ratio of the form fac- 
tors, the fit assumes flat distributions [19]. The result 
of the fit is given in Eq. 3 and the bound on the (p, fj) 
plane is shown in Fig. 3. 



\V td \/\V ta \= 0.16 ±0.02 



(3) 



The most recent value of this ratio obtained from the 
Am s measurement by the CDF Collaboration is about 
two standard deviations away from this result [20]. 
However, it can be argued that the non-penguin con- 
tribution to the charged decay mode B + — » p + j is 
significant. If the fit is repeated only with the neutral 
mode, the new result becomes |Vtd|/|^ts| = 0.16±0.04. 



BABAR results do not agree well with each other for 
the decay B° —> p°j , the clear signal in Belle mea- 
surement of this mode gives the first observation of 
the b — ► e?7 process. Figure 2 shows the beam-energy 
constrained mass M bc = ^ '(££ eam / c 2 ) 2 - |p^/c| 2 , 
and the energy difference AE = E* B - E£ eam distri- 
butions for the reconstructed B-meson candidates. 
Here, pg and E* B are the cm. momentum and energy 
of the B candidate, and El eam is the cm. beam 
energy. The simultaneous fit to all three decay 
modes, assuming isospin symmetry 2 yields the final 



2 It should be noted that the individual fit results are in 
marginal agreement with the isospin relation. From Monte 
Carlo pseudo-experiments assuming the isospin relation, the au- 
thors of the Belle paper find that the probability to observe a 
deviation equal to or larger than the measurement is 4.9%. 



4.6^ sll 

4.1. Inclusive b — > sll 

In addition to the photon penguin diagram of b — > 
57 , at the lowest-order b — > sll process has contribu- 
tions from two other diagrams in the SM: a weak pen- 
guin diagram, where the Z-boson replaces the photon 
and a IF-box diagram. Because of these contributions, 
OPE for b — » sll has ten operators (not counting the 
operators that are suppressed by increasing orders of 
1/mw)- The investigation of b —* sll can identify the 
sign of the Wilson coefficient C7 of the photon penguin 
operator, and can provide the values for the Wilson 
coefficients Cg and C10 of the vector and axial-vector 
electroweak operators. 

The inclusive b — > sll rate has been measured 
by BABAR and Belle, using a sum-of-exlusive-states 
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Figure 3: The bound on the (p, n) plane from the fit to 
the ratio B{B ^ (p/u)j)/B (B — » K* 7) [19]. 



observed exclusive b — > dj decays, 3 they are experi- 
mentally cleaner and it is possible to measure ratios 
and asymmetries in the rates, which can usually be 
predicted more reliably than the branching fractions. 



• Belle 

□ AM et al. 

□ Zhong et al. 



K*i + r 



1 1.5 

Branching Fraction 



2.5 



Figure 4: Comparison of Belle [24] and BABAR measure- 
ments [25] of the B -> Kit and B — ► K* U branching 
fractions with two predictions from the SM, using differ- 
ent form factor calculations [23]. (The plot is courtesy of 
J. Berryhill.) 



(semi-inclusive) method. (So far, unlike b — > sj , 
there has not been any fully-inclusive measurements.) 
The BABAR measurement, based on 82fb _1 of data, 
reconstructs 10 hadronic final states, whereas the 
Belle measurement utilizes 140fb _1 of data and re- 
constructs 18 hadronic final states [21]. The results 
of these measurements, given as a function of the 
dilepton mass squared q 2 , have been compared with 
the next-to-next-to-leading-log predictions for the two 
choices of the sign of C7. As listed in Table III, the 
sign of C7 is quite unlikely to be different from that 
in the SM [22]. 

Table III BABAR and Belle measurements for the inclusive 
branching fraction (in units of 10 ) of B — » X s tt in 
two ranges of the dilepton invariant mass squared. The 
weighted average is compared to the predictions in the 
SM and with reversed sign of C 7 [22]. 



q 2 range 


1 - 6 GeV 2 


> 4ml 


BABAR 
Belle 


1.8 ±0.9 
1.5 ±0.6 


5.6 ±2.0 
4.1 ± 1.1 


Average 


1.6 ±0.5 


4.5 ± 1.0 


SM 

C 7 -» -Or 


1.6 ±0.2 
3.3 ±0.3 


4.4 ±0.7 
8.8 ± 1.0 



4.2. Exclusive b —> s££ 

Two exclusive b — > sil decays have been stud- 
ied extensively at the B Factories: B — ► Kli and 
B — > K* li , where I — e,fi. While their branch- 
ing fractions (Fig. 4) are comparable to the recently 



_l 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 


11 1 j 1 1 1 1 1 j 1 1 i_ 
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» " ;,7 

1 1 1 1 1 1 1 1 


L.,^.,,. n --- y j 
1 1 f 1 1 1 1 1 1 1 1 1 1 1- 



2 4 6 8 10 12 14 16 18 20 



q GeV 2 /c 2 

Figure 5: Forward-backward asymmetry in B — * K* U 
decays as measured by the Belle Collaboration in five bins 
of q 2 [26] . The fit result for negative A7 solution (solid) is 
shown along with some alternative choices for the signs 
and values of the Wilson coefficients; A7 positive case 
(dashed), A10 positive case (dot-dashed), or both A-j and 
j4io positive case. 

Of particular interest is the forward-backward 
asymmetry, Afb, defined for the angle of the posi- 
tively charged lepton with respect to the flight direc- 
tion of the B° or B + meson in the dilepton rest frame. 
In B — > K* tl decay, the SM predicts a distinctive 
pattern for this quantity as a function of q 2 , but the 
presence of new physics can alter its sign and magni- 
tude dramatically. A fit to Apg can be used to extract 
the sign of C7 and the values of Cg and do. 

Figure 5 shows the Afb measurement from the 
Belle Collaboration, in five bins of q 2 . The plot also 
shows the projection of the fit that has been performed 



3 Note that B — * KU is currently the lowest-branching- 
fraction B decay measured at better than 5cr significance. 
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by fixing the value of A7 to its SM value of —0.330 
and choosing Ag/A? and Ai /A 7 as fit parameters, 
where Ai are the leading terms in Cj. Not shown is 
the alternative fit in which A 7 is set to its sign-flipped 
value. The results of both fits, as given in Table IV, 
are consistent with the SM predictions. New physics 
scenarios with positive AgAio are excluded at 98.2% 
confidence. 

Table IV The ratios of the Wilson coefficients obtained 
from the fits to the Belle Afb shape in B — > K* 11 de- 
cays [26]. Two independent fits have been performed, with 
the coefficient A 7 fixed to either its SM value (negative A 7 ) 
or fixed to its sign-flipped value. Also listed are the SM 
predictions for the measured ratios. 





SM 
value 


fitted value 
with A 7 = -0.33 


fitted value 
with A 7 = +0.33 


A 9 /A 7 
A 10 /A 7 


-12.3 
12.8 


-15.3±2;S ± l.l 
10-3±|"J ± 1.8 


-16.3^;? ± 1.4 
11.1±3.9 ± 2.4 



Figure 6 shows the results of the same measurement 
from the BABAR Collaboration. Similar to the Belle 
results, the large asymmetry observed in the high-q 2 
range disfavors new physics scenarios in which the sign 
of the product C9C10 is flipped from its SM value. For 
the low-q 2 range, a lower limit on the asymmetry is 
obtained: A FB > 0.19% at 95% C.L. It should be 
noted that the SM prediction for this range is slightly 
below this lower limit [25] . 





















1 ! 1 1 1 1 ! 1 1 



2 4 6 8 10 12 14 16 18 20 

9 2 (GeV 2 /c 4 ) 

Figure 6: Forward-backward asymmetry in B — > K* ££ 
decays as measured by the BABAR Collaboration in two 
bins of q 2 [25]. The arrow in the low-q 2 bin indicates the 
allowed region at 95% C.L. The SM prediction (solid) is 
shown, as are the predictions from C 7 sign-flipped (dot- 
ted), C9C10 sign-flipped (dashed) and both C 7 and C9C10 
sign-flipped cases. 

In addition to the Afb, BABAR Collaboration has 
also measured the fraction of the longitudinal polar- 
ization, Fl , of K* . Fl is sensitive to left-handed cur- 
rents with complex phases or to right-handed currents 



in the photon penguin amplitude. This measurement 
has the potential for determining the sign of C7 inde- 
pendently, however the current value is in agreement 
with both signs. With 1 ab _1 of data, wrong-sign C7 
could be excluded. 

For B — » Kil , the Afb prediction is zero for all q 2 , 
both in the SM and many of its extensions. Measure- 
ments from both experiments are consistent with this 
prediction. In this decay, BABAR has also attempted 
to extract the fraction of (pseudo-)scalar contribu- 
tion to the decay amplitude, however with the limited 
statistics the result is inconclusive. With more data, 
the developed technique will prove to be useful. 

5. Conclusion 

Radiative penguins have been one of the corner- 
stones of the search for new physics in the last decade. 
While the current measurements have not indicated 
deviations from the Standard Model, the data col- 
lected at the B Factories have significantly improved 
the experimental precision, putting more and more 
stringent constraints on the physics beyond the SM. 
The Wilson coefficients are now experimental observ- 
ables. 

The developed techniques are highly scalable. The 
measurements summarized in this review utilize only 
one-fourth to one-half of the data that is already avail- 
able. Both B Factories continue running with increas- 
ing luminosities and the collected data is likely to be 
doubled before the end of their runs. 

In addition to the direct role they play in the 
search for new physics, radiative penguins also provide 
highly-valueable auxiliary information. The heavy 
quark expansion parameters and the shape function 
derived from inclusive b — > S7 has already been used 
to reduce the uncertainties in the extraction of \V c b\ 
and \V u b\- The measurement of the b — > d'j penguin 
can be considered like the Am s of B Factories. 



Acknowledgments 

I would like to thank the organizers of the con- 
ference for an interesting and enjoyable conference, 
Asoka de Silva for his input on the typesetting of 
this document, and many collegues from the BABAR 
Collaboration for very helpful suggestions during the 
preparation of the talk. 



References 

[1] BABAR Collaboration, B. Aubert et ai, Nucl. lu- 
strum. Methods A 479, 1 (2002); Belle Collabo- 



fpcp06_344 



6 



Flavor Physics and CP Violation Conference, Vancouver, 2006 



ration, A. Abashian et al, Nucl. Instrum. Meth- 
ods A 479, 117 (2002). [16 

[2] T. Hurth, Rev. Mod. Phys. 75, 1159 (2003). 

[3] CLEO Collaboration, R. Ammar et al., 

Phys. Rev. Lett. 71, 674 (1993). [17 

[4] E. Barberio et al., hep-ex/0603003; online update 

at http://www.slac.stanford.edu/xorg/hfag . [18 

[5] Time-dependent CP violation measurements are 
covered in H. Kakuno's talk "New Physics 
Searches at B Factories" at this conference. [19 

[6] Belle Collaboration, K. Abe et al., "Observation 
of B+ -> pAi' , submitted to ICHEP 2004, Bei- 
jing, China, hep-ex/0409009. 

[7] Belle Collaboration, M.-C. Chang et al, 

Phys. Rev. D 71, 072007 (2005). [20 

[8] Belle Collaboration, Phys. Rev. Lett. 94, 111802 

(2005) and references therein. [21 

[9] BABAR Collaboration, B. Aubert et al, "Mea- 
surement of branching fractions and mass spec- 
tra of B — > Kirirj v , contributed to 2005 Inter- 
national Europhysics Conference on High Energy [22 
Physics, Lisboa, Portugal, hep-ex/0507031. 
[10] M. Gronau et al, Phys. Rev. Lett. 88, [23 
051802 (2002); M. Gronau and D. Pirjol, 
Phys. Rev. D 66, 054008 (2002). 
[11] Belle Collaboration, S. Nishida et al., [24 
Phys. Lett. B 610, 23 (2005); BABAR Col- 
laboration, B. Aubert et al, hep-ex/0603054, 
submitted to Phys. Rev. Lett. 
[12] A.J. Buras and M. Misiak, 

Acta. Phys. Polon. B33, 2597 (2002). [25 
[13] T. Hurth et al., Nucl. Phys. B 704, 56 (2005). 
[14] R. Faccini et at, "BABAR Physics Reach", [26 
http://babar-hn.slac. Stanford. edu:5090/hn/ 
aux/ rfaccini/physReach/2dhmodel.html . 
[15] O. Buchmiiller and H. Flachcr, Phys. Rev. D 73, 



073008 (2006). 

E. Barberio, "Inclusive semileptonic B decays", 
contributed to Flavor Physics and CP Violation 
2006 conference, hep-ex/0605098. 
Belle Collaboration, K. Abe et al, hep- 
ex/0506079, submitted to Phys. Rev. Lett. 
BABAR Collaboration, B. Aubert et al, 
Phys. Rev. D 72, 091103 (2005); 
Phys. Rev. Lett. 94, 011801 (2005). 
UTfit Collaboration, M. Bona et al, 
http: / /utfit. romal.infn.it/btovg/ckm-btovg. html. 
This fit makes use of recently improved Light- 
Cone-Sum-Rules calculations of P. Ball and 
R. Zwicky, JHEP 0604, 046 (2006). 
Talk by Guillclmo Gomez-Ceballos, "Update on 
CDF B s Mixing" , at this conference. 
BABAR Collaboration, B. Aubert et al, 
Phys. Rev. Lett. 93, 081802 (2004); Belle Col- 
laboration, M. Iwasaki et al, Phys. Rev. D 72, 
092005 (2005). 

P. Gambino, U. Haisch and M. Misiak, 
Phys. Rev. Lett. 94, 061803 (2005). 
A. Ali et al, Phys. Rev. D 66, 034002 
(2002); M. Zhong, Y. L. Wu, and W. Y. Wang, 
Int. J. Mod. Phys. A 18, 1959 (2003). 
Belle Collaboration, K. Abe et al, "Measurement 
of the Differential q 2 Spectrum and Forward- 
Backward Asymmetry for B — > K^i+i^" , con- 
tributed to Flavor Physics and CP Violation 2004 
conference, Daugu, Korea, hep-ex/0410006. 
BABAR Collaboration, B. Aubert et al, 
Phys. Rev. D 73, 092001 (2006). 
Belle Collaboration, A. Ishikata et al, hep- 
ex/0603018, submitted to Phys. Rev. Lett. 



fpcp06_344 



